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Surface Fermi arcs are the most prominent manifestation of the topological nature of Weyl 
semimetals. In the presence of a static magnetic field oriented perpendicular to the sample sur¬ 
face, their existence leads to unique inter-surface cyclotron orbits. We propose two experiments 
which directly probe the Fermi arcs: a magnetic field dependent non-local DC voltage and sharp 
resonances in the transmission of electromagnetic waves at frequencies controlled by the field. We 
show that these experiments do not rely on quantum mechanical phase coherence, which renders 
them far more robust and experimentally accessible than quantum effects. We also comment on the 
applicability of these ideas to Dirac semimetals. 


I. INTRODUCTION 

Topology in various guises plays a central role in mod¬ 
ern condensed matter physical. In recent years, a sharp¬ 
ened understanding of the topology of electronic wave 
functions in crys tals has stimulated the discovery of new 
phases of matter®!. Among the remarkable manifesta¬ 
tions of these phases are robust gapless edge modes and 
precisely quantized bulk response functions, linked by the 
celebrated bulk-boundary correspondence. Although the 
original applications of topological ideas to band struc¬ 
tures relied on the existence of a fully gapped bulk spec¬ 
trum, more recently it has been recognized that protected 
surface states can arise even in gapless systems-^ 11 . 

The prototypical example of a gapless topological 
phase is a Weyl semimetal (SMjHm 

a three dimensional 
crystalline material where the bulk is gapped except at 
an even number of points in the Brillouin zone (BZ) in 
which the energy bands touch — the Weyl nodes. In the 
vicinity of these nodes, the electrons disperse as mass¬ 
less relativistic particles, and may be described at low 
energies by a Weyl Hamiltonian familiar from particle 
physics, H ~ ±hvok ■ a. Here a is a pseudo-spin degree 
of freedom, k is measured from the Weyl node, and the 
sign is set by the electron chirality 12 . From the perspec¬ 
tive of band topology, a Weyl node is either a source or 
a sink of Berry flux, depending on the chirality. While 
the total Berry flux in the 3D BZ must be zero, as long 
as Weyl nodes of opposite chirality are separated in mo¬ 
mentum space, there must exist two-dimensional cuts of 
the 3D BZ in which the Chern number is non-zero. Each 
such cut defines a 2D integer quantum Hall state. For a 
finite sample, these states necessarily have edge modes on 
the appropriate surfaces. Consequently, for 3D samples 
these real-space surfaces host “Fermi arcs” of states that 
belong to the two dimensional momentum-space Fermi- 
surfac^. 

When a magnetic field is applied perpendicular to real- 
space surfaces that carry Fermi arcs, electrons traverse 
unique cyclotron orbits that connect opposite surfaces of 
the sample. These cyclotron orbits are central to our dis¬ 
cussion. They are reflected in the quantum Shubnikov- 
deHaas (SdH) oscillations of the resistanc^l, as was 


recently observed in a closely related Dirac semimetal 
CdaAsJSl. In addition, Weyl semimetals have been pre¬ 
dicted to exhibit various unusual magnet o-tran sport phe¬ 
nomena, related to the “chiral anomaly T 1 H 1 6 [ 

Following an early suggestion 17 that certain iridium 
pyrochlores may host a semimetal with Nyp = 24 Weyl 
nodes, the number of Weyl SM candidates has prolif¬ 
erated to include HgCr 2 Se 4 with N\y = 'f® and het¬ 
erostructures of normal and magnetically doped topo¬ 
logical insulators 10 . In particular, non-centrosymmetric 
transition-metals, such as TaAs, have been predicted 
in Ref. flUland[2011o be Weyl SMs. Indeed, recent photoe¬ 
mission and transport measurements provide strong evi¬ 
dence for realization of a Weyl SM phase in Ta AfPBEH. i n 
addition, closely related Dirac semimetals, that also carry 
Fermi arcs, have been observed experimentalljf^MSol. 

In this work we show how inter-surface cyclotron or¬ 
bits affect the electronic properties of Weyl semimetals 
already at the semi-classical level. As a result, we are able 
to propose two experiments to probe these trajectories 
without requiring quantum mechanical phase coherence. 
These experiments pose far less stringent requirements 
than SdH oscillations in terms of sample purity, surface 
roughness and temperature. 

We consider a box-shaped slab of Weyl semimetal with 
the main axes being the Cartesian axes, and with the 
Fermi arcs on the z = 0, L surfaces. We assume L to be 
much smaller than the other two dimensions. A magnetic 
field B is to be applied in the z-direction. For concrete¬ 
ness, we consider a time-reversal symmetric Weyl SM. 
Therefore, the surface must include an even number of 
Fermi arcs, and the minimal number of Weyl nodes is 
N\ v = 4, as shown in Fig. EK For simplicity, we ignore 
the curvature of the arcs and consider straight Fermi lines 
directed in the y direction with a constant band velocity 
v a along the x direction and momentum extent k 0 . Al¬ 
though the electrons that reside on a single arc posses a 
non-zero mean velocity, the surface current is zero, due to 
the cancellation between any pair of time reversed arcs. 

We propose two related experiments. In the first, we 
consider two parallel line-shaped Ohmic contacts placed 
on the z = L surface of the slab, separated by a <C L as 
depicted in Fig.[2ji. In standard metals, the current path 
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FIG. 1. (a) Fermi arcs in the surface Brillouin zone for a 
time-reversal-invariant Weyl SM. For simplicity, we ignore the 
curvature of the arcs and consider straight arcs that are di¬ 
rected along the y- axis. The extent of the arcs is ko. The 
red dots denote pair of Weyl nodes with a positive/negative 
chirality. The blue dots denote their time-reversal partners. 
The electron velocity, v a , is perpendicular to the arc at each 
point. Within the straight arcs approximation, the magnitude 
of |v a | = v a is independent of k. In the presence of a perpen¬ 
dicular magnetic field, the surface electrons ’slide’ along the 
Fermi arcs towards the negative chirality Weyl node. An elec¬ 
tric field pulse in the x-direction leads to a momentum shift, 
A k, i.e., the right arc is slightly populated while the left arc 
is depopulated, (b) Bulk LL near one of the Weyl nodes. The 
red line denotes the chiral LL. 


extends a distance of the order of a into the bulk. Hence, 
for a -C L the current at the z = 0 surface vanishes as 
a/L 2 . We show that in the presence of a perpendicular 
magnetic field when a voltage V is applied between the 
two contacts on the top surface, opposite currents are 
induced in the two surfaces of the sample. As a conse¬ 
quence, a voltage aV is induced in the bottom surface, 
that scales linearly, a ex |B|, for small fields. 

Our second proposal considers an electromagnetic mi¬ 
crowave radiation propagating from z = oo downward 
along the - 2 -axis, as depicted in Fig. §*■ We show that 
when the slab is much thicker than the semimetal skin 
depth, such that the radiation is expected to be mostly 
reflected by the slab, there are transmission resonances 
at which a significant part of the radiation is transmit¬ 
ted, with an amplitude that is independent of L. The 
transmission amplitude is again linear in |B| for small 
fields, as is the resonant frequency. 

The effects we discuss here involve transfer of elec¬ 
trons between Weyl nodes. In a clean Weyl semi-metal 
in a magnetic field, there are two mechanisms for elec¬ 
trons to be transferred between different nodes: the chiral 
anomaly (which is effective in the bulk) and the Fermi 
arcs at the surface. Our proposals rely on the second 
mechanism, in contrast to those of Ref. [TB] that originate 
from the first. 


FIG. 2. (a) DC configuration: a current, I, is injected via 
the upper left contact and withdrawn from upper right con¬ 
tact. The voltage difference between the lower contacts is 
measured. Without the non-local orbits, current flow in a 
’cigar’ shape configuration between the contact, (b) AC con¬ 
figuration: current cycle due to a pulse of EM field on the top 
surface leads to emission from the lower surface. 

II. NON-LOCAL CONDUCTIVITY 

Both these phenomena originate from the same mech¬ 
anism — the non-local conductivity of a Weyl semimetal 
in a magnetic field. Linear response theory defines the 
conductivity in space-time through the relation j(r. t) = 
J dr'dt'Y,(r, r', t — F)E(r , ,f'). Applying Drude theory to 
a doped Weyl/Dirac node, we find 

E = Ed = ^-uofc)U -t/T 6(r - r') = a 0 S( r - r'), (1) 

where r is the momentum relaxation time, vq the velocity, 
kF the Fermi momentum and do is the local Drude con¬ 
ductivity. We use E to denote time-domain conductivity, 
and remove the tilde for frequency-domain conductivity. 

The unique inter-surface cyclotron orbits result in a 
non-local contribution to the conductivity E, not cap¬ 
tured by Eq. 0. Consider an electric field pulse induced 
by a vector potential A x (z, t) = —A 0 f(z)6(t), where f(z) 
is spread over a short length scale d around z = L. This 
length scale is a few times the spatial extent of the surface 
states in the 2 -direction. Immediately after the pulse, a 
surface current emerges, as may be computed from the 
shift of the Fermi level of the two Fermi arcs due to the 
applied vector potential. The induced current density is 
jo = ev a koAk where HAk = eAof(z). 

Following this shift of the Fermi surface, the magnetic 
field B causes the surface electrons to ‘slide’ along the 
Fermi arcs towards the negative chirality Weyl point, at 
a rate k = |v a x B (Fig. 1^,). As electrons slide on 
the Fermi arcs they eventually arrive at the Weyl nodes, 
where they merge into the 3D bulk. In the presence of 
the magnetic field, the bulk spectrum in the vicinity of 
each Weyl node consists of dispersive Landau levels (LL) 
as depicted in Fig. [T}>. Specifically, the chiral LL that 
originates from the v = 0 LL of Dirac electrons serves 
as a direct passageway for electrons from the top to the 
bottom surface. Note that the electrons travel along the 
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same cyclotron orbit responsible for the unique quantum 
oscillation signature of Weyl SMs 13 . 

When the excess current is all in the chiral state , i.e., 
when the Fermi energy p satisfies \p\ < voV2HeB, the 
entire current flows to the opposite surface, where it flows 
in the — x direction. In the z = 0 layer the sign of the 
velocity reverses, and hence so does the sign of k y , such 
that the electron motion along the arc brings it to a chiral 
state that flows back to z = L, completing the inter-layer 
cyclotron motion. In the absence of scattering, this cycle 
repeats indefinitely. Scattering between different Weyl 
nodes is detrimental to inter-surface cyclotron orbits. In 
clean samples such scattering is rare due to the large 
momentum difference between nodes. Scattering within 
a Weyl node is highly suppressed when |^/| < voV^heB 
since there are no Landau levels to scatter into; we discuss 
the case |p| > v 0 V2heB below. 

Altogether, then, a pulse of an electric field in the in¬ 
direction at z = L leads to alternating and opposite cur¬ 
rents j \, j 3 in the two surfaces, and an alternating bulk 
current j -2 — ji between the surfaces. The periodic or¬ 
bit is depicted in Fig. The period To is obtained by 
combining two basic timescales: T 0 = 2(T z +T arc ), where 
T z = L/v o is the time needed for electrons to cross from 
the upper to the lower surface, and T arc = hko/{eBva) 
is the time needed for electrons to slide along the Fermi 
arc. Here again, kg is the extent of the arc in /c-space 
and v a is the magnitude of the arc velocity as depicted 
in Fig. [l}i. 

We now focus on the current in the lower surface, j '3 in 
Fig. [ 2 ^, from which we calculate E(0, T, t — t'). We con¬ 
sider only the evolution of the excess electrons due to the 
imbalance generated by the electromagnetic (EM) pulse, 
since the current is solely determined by these electrons. 
At t = 0 the excess electrons populate all the states along 
the Fermi arc. These states are gradually depleted to the 
bulk at a constant rate, j^Bv a , until all states are empty 
at t = T arc . The first excess electron reaches the lower 
surface at t = T z . Gradually, more and more excess elec¬ 
trons reach the lower surface until t = T z + T arc , where 
a maximum in the current density, j 3 , occurs. Then, the 
excess electrons start to leave the lower surface and move 
towards the upper one, completing the cycle. Scattering 
between the Weyl nodes is expected to suppress this cur¬ 
rent, and we characterize it by a relaxation time t v and 
length l v . The current decays as e~ t ^ Tv as more and more 
cycles occur. The effect of intra-node scattering, which 
takes place when |/i| > v 0 \j2heB is discussed below. 

A sketch of the real-time current on the lower surface, 
j 3 (f), is shown in Fig. [ 3 ^, for the case \p\ < v 0 \/2HeB. 
Assuming a linearly rising current as the electrons arrive 
at the bottom surface, the current on the lower surface 
is obtained by convolving a periodic function represent¬ 
ing the cycles, with a ‘triangle function’ A (t) describing 
the growth and decay of the current in each cycle, and 



FIG. 3. (a) A sketch of real time current on the lower surface 
for |/r| < voy/2heB. In this case, the source of current decay 
is inter-valley scattering, (b) The transmission coefficient as 
a function of the source frequency. A significant part of the 
radiation is transmitted whenever the frequency of the applied 
EM field is an integer multiple of u>o = 2 tv/To- 


scattering induced exponential suppression. 

h{t) ~ A(t) * ^2 6 (t-T 0 (n+ e~^ (2) 

where >}= denotes a convolution. The current induced at 
the same surface at which the electric field is applied is 
given by a similar expression, with the 1/2 absent. 

Fourier transforming and substituting jo we get an ex¬ 
pression for contribution of the cyclotron orbits to the 
conductance, 


£(0 ,L,u>) 
£ (L,L,u) 


-jx{ 0,tu) 

E x (L,u>) 

jx{L, u}) 

E x (L,u) 


iG 0 sinc 2 ^ uT * IC ^ 
2 sinh {ip) 

E(L,0,w)e* 


( 3 ) 

( 4 ) 


where <p = (t v 1 — ico)To/2 and Go = e 2 u a fcoT arc //id. Con¬ 
sequently, we approximate the non-local currents by 


j^ L « d*w(u)6(z - L ) [E x {L)e* - E x { 0)] + (0o L), 

( 5 ) 

where a w = £(0 ,L,u) is the ‘Weyl conductivity’. The 
Drude conductivity <to and the Weyl conductivity aw 
have the same units. In our estimate of aw we neglect 
the variation along the arc of the localization length of 
surface states in the ^-direction. We also ignored diabatic 
transitions into the bullP^. The latter effect is just a 
correction fco —> ko — fM^ 1 where £b is the magnetic 
length and (3 ~ 0(1). 


III. EXPERIMENTAL SIGNATURES OF 
NON-LOCAL CONDUCTIVITY 

Having calculated the non-local part of the conductiv¬ 
ity, we are in a position to analyze the two experiments 
we propose. 
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A. DC Transport 

When a DC voltage is applied between the lines x = 
±a/2 on the z = L surface, a current flows in the sam¬ 
ple. The current density and electric field must satisfy 
Kirchoff’s rules and Ohm’s law, 


Vj = 0 

(6) 

V x E = 0 

(7) 

j = £E 

(8) 


where £ is the calculated conductivity (including the 
non-local part in Q). We ignore the inter-valley cur¬ 
rents due to the chiral anomaly, that would generically 
make the bulk conductivity moderately anisotropic (see 
the next section for a discussion). The boundary condi¬ 
tions impose zero current perpendicular to the surface ev¬ 
erywhere except at the contacts, and enforce the voltage 
V = - / E • dl between the contacts. At zero magnetic 
field the conductivity is purely local. Consequently, the 
current path extends a distance of the order of a into the 
bulk, and vanishes as a/L 2 at the z = 0 surface. Fig. [ 4 ] 
presents the current path in the presence of a non-local 
conductance. The color represents the stream function 
ip(x, z) which is related to the current as follows: 



x x 


(a) (b) 

FIG. 4. The stream function ip for (a) pw/po = 0 and (b) 
Pw/Po = ovc/c 0 = 0.15. The current flows along contours 
(green) of constant ip (in the direction of the arrows), and 
its magnitude is given by | Vgb| which is proportional to the 
density of the green contours. In (a) we find that ~ 0.005 
while in (b) we find Wa « —0.273. 


This is the so-called chiral anomaly. The change in the 
bulk current is then given by: = evo^. Employing 

theses relations, we conclude that an extra current j z is 
produced due to the chiral anomaly: 



Jx = ~d z ip, J z = d x ip. (9) 

A few equal value contours of ip(x, z) and the direction 
of the current along them are also shown in to Fig. [4] 
The current flows along contours (green) of constant ip 
(in the direction of the arrows). The values of the pre¬ 
sented contours are equally spaced. Hence, the current 
that flows between any two contours is the same, and 
the magnitude of the current density is proportional to 
the density of the (green) contours. As can be seen in 
Fig. gji, only a small current reaches the lower surface. 
On the other hand, as seen in Fig. [4 Jd, the application of 
a magnetic field leads to an opposite current at the z = 0 
surface. The chiral modes in the bulk, which mediate the 
current between the two surfaces, are always parallel to 
the magnetic field. Therefore, when the magnetic field is 
along the z-axis, the current on the lower surface appears 
approximately below regions on the upper surface where 
a non-negligible electric field is developed. Details of the 
calculation are given in Appendix | A 1| 


where £3 is the magnetic length. Evidently, the chiral 
anomaly renormalizes the longitudinal conductivity in 

2 r 

the 2 direction, i.e., a zz = ao + while the other 

B 2 O 

directions remain unchanged, a xx = <jq = jpkp£. Hence, 
the effect of the chiral anomaly is to make the bulk con¬ 
ductivity anisotropic. The solution of the Kirchoff-Ohm’s 
equations in the anisotropic case may be related to the so¬ 
lution in the isotropic case by the following simple rescal¬ 
ing: 

(E X ,E Z ) = ( E x ,oT l E z ), 

{.jxijz) = (jx,ajz), 

( x,z) = ( x,az ), and aw = adw 

where the tilde represent quantities evaluated in the 
isotropic case, and a = ( a zz /a xx ) 1 / 2 ~ ((2A + 

1 ) -1 L/£) 1 / 2 with N being the filling factor of the non- 
chiral Landau levels. Thus, as long as a ~ 0(1) the DC 
effect is only moderately affected by the chiral anomaly. 


1. Inter-valley currents due to chiral anomaly 


B. AC Transport 


We address the effect of the chiral anomaly on the DC 
effect. When a DC voltage is applied between the two 
upper contacts, the electric field lines bend into the bulk 
at the vicinity of the upper surface. The resulting E z 
leads to a density transfer between the two chiral modes 
in the bulk according to: 


dn 

dt 



E B = 



E Z B 


The second experiment examines the transmission of 
electro-magnetic radiation through the slab. Assuming 
translational invariance in the xy plane, the propagation 
of a monochromatic electromagnetic wave along the z- 
direction is described by the inhomogeneous wave equa¬ 
tion (c = 1), 


( 10 ) 


(d 2 z + w 2 ) E = - (iw/xj - ( d z p)z ) 


( 11 ) 
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where p and j are the charge and current densities inside 
the metal and p is the permeability of the metal. Assum¬ 
ing incident radiation polarized along x and with an am¬ 
plitude Eq, only the ^-component of the electric field is 
nonzero outside the slab, so that the solution of Eq. (11) 
satisfies E(w, z = 0) = E( 0)x and E(w, z = L) = E(L)x. 
Solving for E( 0) as a function of E 0 and ui allows us 
to characterize the transmission amplitude of the wave 
through the sample, |E( 0 ,w)/E 0 |. 

In the presence of a non-local electrodynamic response, 
we can rewrite Eq. (fTTl) as 


(d 2 + ui 2 + icupcro) E x (z, u>) = i upj 


NL 


( 12 ) 


■NL 


where do is the conductivity of the metal, and j 1-iJ is 


defined in Eq. ([5]). For E = 0, the solution of Eq. (12) 
is an evanescent wave decaying exponentially into the 
bulk, E x (z) ~ e~ z ! & with the skin depth <5 _1 = Im(w 2 + 
iuipao) 1 ^ 2 ~ (cLi/rcro/2) 1 / 2 . We assume that the bulk is 
sufficiently metallic that the Hall angle is small; then, 
in the absence of non-local terms, the transmission is 
exponentially small in the sample width. 

Solving Eq. (12), we find a significant increase in the 


transmission coefficient when the frequency of the applied 
EM field is an integer multiple of to o = 2tt/T 0 , as depicted 
in Fig. !>■ In appendix [Bj we analyze the transmission 
coefficient as a function of the ratio aw /cr 0 . We find that 
for frequencies where a w(w) is non-negligible, electric 
held E( 0) is comparable to E(L). As aw increases, the 
transmission increases rapidly until aw ~ <7o, where it 
saturates. The saturation value depends on the ratio 
ui/ao . The non-local orbits “connect the surfaces”, which 
eventually leads to a held in the lower surface E( 0) that 
may be as high as the held E(L). That does not amount 
to a full transmission because E(L) is not E 0 . For a 
thick standard metal the radiation is partially reflected 
and partially absorbed. The absorption originates from 
the penetration of an electric held into the metal. In the 
Weyl case, the mechanism we described transforms part 
of the absorption into transmission, with an effectiveness 
that increases with aw/&o- However, this mechanism 
does not completely eliminate the rehection, and hence 
the transmission saturates at a value smaller than one. 
The increase in absorption with u> in standard metals is 
converted into resonant transmission at multiples of the 
cyclotron frequency. 

In other words, in the presence of a magnetic held the 
Weyl semimetal acts as a “narrow band-pass filter” for 
incident EM radiation around harmonics of the frequency 
cjq = 27 t/Tq. The frequency ujq can be tuned easily, since 
it depends on both the system width, L, and the mag¬ 
netic held B. In crystals where the surface hosting Fermi 
arcs lacks C 4 symmetry, the results of the AC experiment 
would be polarization dependent. More details appear in 
appendix [B] 

The resonant transmission we hnd resembles the 
Azbel’-Kaner cyclotron resonanc^ 1 ®!!, in which mag¬ 
netic helds parallel to the sample surface lead to enhanced 
transmission when the cyclotron radius and frequency 


match the sample thickness and the frequency of the ap¬ 
plied electromagnetic held, respectively. Note that previ¬ 
ously proposed nonlocal effects in Weyl SMs 1 ^ are quite 
distinct from those considered here. 


IV. EFFECTS OF SCATTERING 

Shubnikov-deHaas oscillations of the resistance are a 
consequence of quantum interference, and are hence sup¬ 
pressed by elastic and inelastic scattering, even at small 
momentum transfer. As such, they are suppressed expo¬ 
nentially in the ratio of the size of the cyclotron orbit to 
the mean free path and the ratio of the temperature to 
the cyclotron energjEH. For the inter-surface cyclotron 
orbits we consider this requirement restricts the sample 
size L < l, achieved through low temperature measure¬ 
ments on very thin, high-purity samples. 

In contrast, the two transport phenomena that we 
discuss here, while they rely on the existence of the 
cyclotron orbits, do not rely on quantum interference 
around the orbits, and therefore phase coherence is unim¬ 
portant. Their existence is affected instead by how mul¬ 
tiple scattering events in the bulk alter the structure of 
E(t) (shown in Fig. [ 3 ^,). To address this, we return to the 
spectrum depicted in Fig. [TJd, and assume the chemical 
potential p to be such that the chiral state overlaps in 
energy with N > 1 Landau levels. 

For concreteness, we consider a Weyl node where the 
chiral state flows downwards (in the — z direction). In 
the absence of scattering the Weyl node has N +1 modes 
flowing downwards and N modes flowing upwards. We 
model the semiclassical scattering as a diffusion process, 
in which the electron is scattered between the 2 N + 1 
modes, with a scattering event taking place at a rate 
1/r, with t «Ij. For each Landau level, the upwards 
and downwards velocities are identical in magnitude, but 
velocities are not necessarily identical between Landau 
levels. Due to the unpaired chiral state, the diffusion 
process results in a downwards average drift velocity Vd, 
which is smaller than the velocity of the chiral mode, Vq- 
Neglecting velocity differences between different Landau 
levels, Vd = vq/(2 N + 1), and the average crossing time 
between surfaces becomes (T z ) = T. > A Moreover, the 
diffusion leads to fluctuations in the crossing time, char¬ 
acterized by a variance AT Z around (T z ). The relative 
width of the distribution of arrival times to the bottom 
surface is given by 

AT Z _^D(X)_ /(2N + l)l 

(13) 

where D = VqT is the diffusion constant, and l = vqt is 
the mean free path. Thus, as long as l(2N + 1 )<L, the 
first passage of the current from one surface to another 
occurs at a rather well defined time. In terms of E (t), 
this will lead to a shift, smoothening and broadening of 
the first peak of Fig. The area under the peak stays 
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constant, since the entire current crosses from one surface 
to another (see below). As the current continues along 
the inter-surface cyclotron orbits the spread of the cross¬ 
ing times increases, and hence the peaks of E(t) further 
broaden, until effectively merging together at a time t 
where Dt ~ L 2 . 

Interestingly, intra-Weyl node scattering does not sup¬ 
press the cyclotron current that gives rise to the non-local 
E(£), even after the peaks merge together. The reason 
for that may be understood by considering an analog sit¬ 
uation in the realm of the two dimensional quantum Hall 
effect. Imagine a v = 1 QHE state at the half-plane 
x < 0, with an downwards-moving chiral edge along the 
y-axis. Now imagine coupling the y < 0 part of the edge 
to N semi-infinite quantum wires of spinless electrons, 
each carrying one upwards-moving chiral mode and one 
downwards-moving chiral mode. Independent of whether 
the coupling is ordered or random, it cannot block a 
downwards moving current emanating from y = +oo, 
due to the chirality of the QHE edge. The coupling to 
the wires merely renormalizes the velocity of the edge 
state and modifies its wave function. 

In Weyl semimetals the surface plays the role of the 
QHE y > 0 edge. As long as the surface is free of scat¬ 
tering, the momentum of states on the surface is well 
defined. The finiteness of the arc, together with the drift 
imposed by the equation of motion hk = ev a x B, enforce 
a flow of the electrons into the bulk, which plays the role 
of the y < 0 region in the QHE case. This flow cannot 
be reversed as long as electrons are not transferred to a 
Weyl node with an opposite chirality, since the direction 
of k is fixed for every arc. Thus, the entire current that 
enters the bulk from one surface must cross all the way to 
the other surface. As long as electrons’ chirality changes 
only at the surface, the peaks in Fig. preserve their 
area under scattering. 

Transforming these observations into the frequency do¬ 
main, we find that E(w) shows clear resonances, as in 
Fig. [3jo, around harmonics of a modified cyclotron fre¬ 
quency 2tt/(T 0 ). The height of the peaks decays fast 
with frequency, but the low frequency peaks remain al¬ 
most unaltered when compared to the clean case. The 
DC part, E (oj = 0), which involves the integral over all 
times, is not altered by intra-node scattering. Conse¬ 
quently, intra-node scattering does not affect the signal 
in the first experiment we proposed here. Two factors 
that do affect this experiment are inter-node scattering 
and the chiral anomaly. Inter-node scattering is charac¬ 
terized by valley relaxation length, l v , which can be tens 
of microntP®. It should therefore be rather ineffective. 
For the AC experiment, inter-valley scattering suppresses 
the area under the peaks in E (£) and hence also the res¬ 
onances in E(w). However, as long as the rate is smaller 
relative to the resonance frequency, its effect is small. 

The effect of temperature on the two experiments is 
indirect, through its effect on the intra- and inter-node 
scattering rates, as well as on the number of bulk filled 
Landau levels. In all cases, the resulting effect scales 


like a power-law. This is in contrast to quantum os¬ 
cillations, whose amplitude decreases exponentially with 
temperatur^H. 

Finally, we comment on the applicability of these ideas 
to Dirac semimetals. Dirac SMs host similar inter-surface 
cyclotron orbits. In Dirac SMs, unlike in Weyl SMs, 
the counter propagating chiral LL reside at the same 
node. This reduces the robustness of the proposed effects, 
since now scatterings between the two chiral channels do 
not require large momentum transfer. Furthermore, in 
Cd 3 As 2 , the surface on which arcs are seen breaks the 
protecting symmetrjJ 13 ! and hence, one might worry that 
the arcs could be reconstructed into an ordinary Fermi 
surface. Nonetheless, as long as the scattering rate be¬ 
tween the different chiral levels is small, as suggested 
irP, the effects should be visible also in Dirac SM. Note 
that the experimental observation of Fermi-arc orbits in 
CdaAsJ-SI i s very encouraging in this regard. 


V. ESTIMATES OF SCALES 


We provide quantitative estimates of various relevant 
quantities. For typical values of L = 5 — 50 /.mi, vq = 
1 • 10 5 — 6 • 10 5 m/s and T arc < T z , the Weyl cyclotron fre¬ 
quency is ui c = 5 — 500 GHz., i.e., in the microwave range. 
In order for the described resonances to be a relevant 
probe, the sample should be thicker than the skin depth 
S at moderate frequencies (if not, samples will be trans¬ 
parent at these frequencies independently of the magnetic 
field.) Assuming low-temperature scattering-rates simi¬ 
lar to the ones measured in Dirac SM, we may e xtract 
a typical sample resistivity of the order ~ 10 /iflarPSE? 7 . 
This value produces a skin depth 5 ~ 40a; -1 / 2 cm, which 
is ~ 1 /im at microwave frequencies, serving as a lower 
bound on sample thickness. The upper bound is deter¬ 
mined by the inter-node scattering processes: we require 
that L < l v , the valley relaxation length, which can 
be tens of micrond 1 ^. As a final point, we remark on 
a subtlety: unlike usual cyclotron orbits, the frequency 
for electron motion on the nonlocal orbit connecting op¬ 
posite surfaces is thickness-dependent. Since the Weyl 
cyclotron resonances occur at the harmonics of this fre¬ 
quency, in order to observe the resonant transmission, 
we require that the skin depth at the lowest resonant fre¬ 
quency is smaller than the sample size, L, which induces 
an implicit dependence of the skin depth on the thick¬ 
ness, S = 5{u c {L)\ in order to observe Weyl cyclotron 
resonance. Thus, we have an additional constraint on 
the thickness, 5[uj c (L)\ < L. Fig. [5] shows the depen¬ 
dence of S on the thickness for various values of the drift 
velocity Vd\ we see that for L > 5/xm, this condition is 
indeed satisfied. 
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FIG. 5. Dependence of skin depth at the lowest resonance 
frequency, 8 \ui c {L)\, on sample thickness L and for po = 
10/rflcm; for L > 5 pm, 8 <C L. 


VI. CONCLUSIONS 

We have shown that inter-surface cyclotron orbits af¬ 
fect the electronic properties of Weyl semimetals already 
at the semi-classical level, and we proposed two exper¬ 
iments which directly probe these orbits. We demon¬ 
strated that, in the presence of a magnetic field, the 
transport in Weyl SMs is characterized by a non local 
conductivity which leads to the appearance of a mag¬ 
netic field dependent non-local DC voltage and to a res¬ 
onant transmission of electromagnetic radiation through 
macroscopic samples of Weyl semimetals. Furthermore, 
we argued that the semi-classical origin of the effects 
reduces dramatically the experimental requirements of 
thin-films, low temperatures and pure samples. 
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Appendix A: DC Setup 
1. Solution of the Kirchoff-Ohm’s equations 

We consider a sample with four stripe-like Ohmic con¬ 
tacts lying along the y direction and at a distance a -C L 
from one another as shown in Fig. §i. Two contacts are 
positioned on the 2 = 0 plane and two on the z = L 
plane. We assume that the y direction extent is much 
larger than L, such that we expect translation invariance 
along y. Therefore, we consider a slice of the x — z plane, 
Fig. [6jo. A current, I, is injected via contact A and with¬ 
drawn from contact B. The voltage difference Vqd will 
then be calculated. 


4 T 




FIG. 6. (a) The DC setup, y direction extent is much larger 
than L. (b) x — z slice of the setup in (a). A current, I, is 
injected via contact A and withdrawn from contact B. 

The equations governing the current flow are: 


<1 

Ch 

II 

0 

(Al) 

V x E = 0 

(A2) 

E = E -1 J, 

(A3) 


We choose to work here with the resistivity and not 
the conductivity since it simplifies the implementation 
of boundary conditions in this setup. The generalized 
Ohm’s law can be expressed in a matrix form: J(z) = 
f dz'E(z, z')E(z'), where T,(z, z') = cto< 5 2jZ ' + aw(z, z'). 
Hence, the inverse relation is given by: 

E(z) = T,~ 1 {z,z')J{z') = {&o8 zz > + a w (z,z')}~ 1 J(z') 

» U(«) - ^ A J( z') (A4) 

00 00 

In the last equality we assumed that aw -C 0o- This ap¬ 
proximation is not essential. However, it simplifies the 
calculation. Later, we also discuss the opposite limit 
aw 0o- Using the form of the non-local conductiv¬ 
ity, Eq. |5]), we conclude that: 

E L (Y) = — J (z) = p 0 J(z) 

00 

E nl (z) « —dp w S(z - L) [J(0) - J(L)] +(0o L), 

















where the L and NL refer to local and non-local, re¬ 
spectively, and pw = &w ( w = 0)/crg. Notice that 
Pw/Po = a w /a o- 

Introducing the stream function ip{x, z) as follows: 


J = y x Slip -t J x = ~d z ip, J z = d x ip. (A5) 


Eq. (Al I and the x , z components of Eq. (A2) are auto¬ 
matically satisfied. The y components of Eq. (A2): 


0 = (V x E) y = (V x E L ) y + (V x E NL ) y (A6) 

= p 0 VV + (V x E NL ) y p 0 V 2 V> = -(V x E NL ). y 


Since no current can leave or enter the system except 
at the contacts, the stream function ip must be con¬ 
stant along the boundaries. Since a current I is injected 
at contact A and withdrawn at contact B, integrating 
J x = —d z ip across the step discontinuity at either con¬ 
tact implies that the jump in ip across contact A ( B ) is 
I{—I). We choose 7=1, thus, the boundary conditions 
are: ip = 1 on the boundary segment between A and B, 
and ip = 0 on the all other boundary segments. Notice 
that the current density J lies along contours of constant 
ip, with magnitude \Vip\- Additionally, the voltage dif¬ 
ference between any two points on the boundary is given 
by: 


Vij — 




Xj 

- J dxp 0 d z ip (A7 ) 

Xi 


For E nl = 0 the solution is depicted in Fig. [ 4 ^ 1 . Evi¬ 
dently, the current path extends a distance of the order 
of a into the bulk, and is therefore vanishingly small at 
the z = 0 surface. Inserting the solution in Eq. (A7) 
gives: V C d/Vab ~ 0.005. 

Next, we introduce the non local part E NL . The 
boundary conditions remain unchanged. The solution 
to Eq. (|A6| with pw/po = aw/a 0 = 0.15 is depicted 
in Fig. gf>, and the voltage ratio becomes: Vcd/Vab ~ 
—0.273. 


In the limit aw(u = 0) <C ao, in which the voltage 
ratio may be calculated numerically, it is proportional to 
aw/a o- The voltage ratio is depicted in Fig. [ 7 } 

Next, we comment on the limit % cr 0 . Neglect¬ 
ing the dependence of the current on x, we may write a 
one-dinrensional model for the current in the ^-direction. 
By discretizing the z-coordinate 2 = (0, A z, 2 A z, ■■ ■ ,L), 
where A z is larger than the thickness of the Fermi arcs, 
the conductivity may be expressed as a finite matrix of 
the form: 


"<To -h aw 0 0 — aw 

0 do 0 0 

0 0 Co 0 

. ~&w 0 • ■ • 0 a 0 + &w- 


(A8) 


Inverting the conductivity matrix yields the resistivity 
matrix: 


-A 0 • • • 0 Br 

0 1/(T 0 0 0 


0 0 1/cto 0 

B 0 • • • 0 A. 


(A9) 


Where A =2k + 2(a 0 +2, w ) and B = *o(Z+2* w) - 


In agreement with Eq. (A4), for aw "C op the co¬ 
efficients become A ss l/a 0 and B « owl<x§- I n 
the opposite limit, aw <r 0 , the coefficients become 

A k, B k, l/(2cro). Independently of the current density 
profile, the condition A ss B forces the electric fields on 
the top and bottom surfaces to be equal, i.e., the voltage 
ratio tends to unity in the limit aw ^ &o- 


Finally, we comment on the locality of E in the x — y 
plane. In the above calculations, we assumed that E is 
local in the x — y plane. The locality breaks for length- 
scales smaller than fc 0 1\. Nonetheless, as long as kol 2 B <C 
a the effects of the non-locality in the x — y plane are 
negligible and the locality assumption is justified. 


Vcd 

Vab 



0.1 0.15 

aw(uj=0) 

(To 


0.25 


FIG. 7. The ratio of the voltages on the bottom and top 
surfaces as a function of the ratio aw/ao. The voltage ratio 
is calculated numerically, in the limit aw/ao ^ 1, using the 
approximated Eqs. |Al||A3| In the limit aw ^ ao, the voltage 
ratio tends to unity. 
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&w /&0 


FIG. 8. The transmission coefficient as a function of S = 
aw /do for different values of fi = (fch) -1 and for A = L/5 = 
20 . 


Appendix B: AC setup — Calculation of the 
transmission coefficient 

The Maxwell equations in SI units are: 


„ „ <9H . TT 

V x E = —ii—— = iwuH 

at 

<9E 

V x H = J + e-^— = J — iweE 

at 


(Bl) 


where /i = iio/i r and e = e 0 e r are the permeability 
and permittivity of the metal. We choose the applied 
a.c. electric field to be along the x directi on. Hence, we 
may conside r on ly the x component of Eq. (Bl). The cur¬ 
rent in Eq. (Bl I has two parts: 3(z) = cr 0 E(z) + J NL (z), 
where do is the conductivity, and J NL (,z) is non-local part 
(see Eq. 5 in the main text). In general, the do depends 
on u>. However, for microwave frequencies cot <C 1, hence, 
we neglect the frequency dependence of d 0 . 


Inserting the two parts of J to Eq. (Bl), assuming 


translational invariance in the xy plane, and intro duc- 


ing the dielectric function, e(w) = e r + Eq. (|Bl|) 
becomes, 




E(z) = iw/hP L (z), 


(B2) 


where we assumed y r ~ 1 and employed the relation 
/i 0 e 0 = c -2 . Integrating Eq. (B2) form 0“ to 0 + and from 


L to L + we find the boundary conditions (BC): E x (z) 
is a continuous function, in particular ^(O - ) = -E x (0 + ), 
E X (L~) = E X (L+), and 


(&zEx)z= o+ (dzEx) z= o - — i ujfid&w [-£+(0) E x (E)\ , 
( d z E x ) z=L + - (d z E x ) z=L - = iuiydaw [E x { 0) - E X (L)]. 

The last two BC are equivalent to the requirement 
that the discontinuity of the tangential magnetic field 
is proportional to the surface current. Next, we consider 
an incoming monochromatic field with a unit amplitude 
Eq = 1. For z ^ 0 or i, the solution to Eq. (B2) has the 
following general form: 


f e~ lkz + E r e ikz , L < z 
E x {z) = i , 0 < z < L , (B3) 

{E t e- ikz , z < L 


where k = ^ and q = ^ In typical metals and for 

microwave frequencies do eo^. In this regime, all the 
physical quantities may be recast in terms of following 
three dimensionless quantities: fl = (kd ) -1 , A = L/5 
and S = aw/co- Imposing the BC, we find 


\Ed 


- 1 - 

E 2 


1-2SQ 

E r 


0 

E t . 


2 SQ 


(B4) 


where the matrix M is given by 

/i i -i o \ 

/ (l+i)n -(l+i)0 1+250 -2SOe ,A/!2 \ 

e (i-DA e -ti-DA o _ e iA/fi 

\-(l+i)Oe (i_1)A (l+i)Oe _(i_1)A -2SO (2SO+l)e iA/n / 


Eq. (B4) may be solved for \E t \ to yield the transmission 
coefficient as depicted in Fig. [8] 

For S = 0, the field in the metal decays as E s e~ z ^ s , 
where 5 ~ (wdo)~ 1,/2 is the skin-depth and the quantity 
E s is the field on the surface. The ratio between the 
field on the surface and amplitude of the applied field is 
proportional to (wdo) 1 / 2 . In particular, the dissipation 
is f dz E • J ~ <5do-E 2 a ( w/d 0 ) 1 ^ 2 . For a fixed do, as 
to increases, the decay becomes stronger, however, the 
value from which it decays becomes larger. 

In the Weyl case, the mechanism we described trans¬ 
forms part of the absorption into transmission, with an 
effectiveness that increases with d^/do- However, this 
mechanism does not eliminate the reflection, and hence 
the saturation of the transmission at a value that is 
smaller than one. The increase of absorption with ui in 
standard metals is transformed here into higher trans¬ 
mission. 
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